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Magnetization measurements were performed to investigate the critical behavior of the 
field-induced magnetic ordering in gapped spin system TICUCI3. The critical density of the 
magnons was obtained as a function of temperature and the magnon-magnon interaction con- 
stant was evaluated. The experimental phase boundary for T < 5 K agrees almost perfectly 
with the magnon BEC theory based on the Hartree-Fock approximation with realistic disper- 
sion relations. The phase boundary can be described by the power law [H N (T) - H c ] oc T 4, . 
With decreasing fitting temperature range, the critical exponent <f> decreases and converges at 
0bbc = 3/2 predicted by the magnon BEC theory. 
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Quantum spin system composed of antiferromagnetic 
spin dimer often shows a gapped singlet ground state. 
In an external magnetic field exceeding the energy gap 
A, S z = 1 component of the spin triplet is created in 
the system. The field-induced S z = 1 component has the 
characteristics of boson and is called magnon or triplon. 
Magnons move to neighboring dimers and interact with 
one another due to the transverse and longitudinal com- 
ponents of the interdimer exchange interactions, respec- 
tively. Consequently, the spin dimer system in the mag- 
netic field can be represented as an interacting boson 
system. 1 Magnons can undergo Bose-Einstein condensa- 
tion (BEC) in a magnetic field higher than the critical 
field H c = A/gfiB, which leads to field- induced magnetic 
ordering (FIMO) . 2 ' 3 Nikuni et al. 2 discussed the FIMO 
observed in TICUCI3, 4 applying the Hartree-Fock (HF) 
approximation to a simplified model 
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where is the kinetic energy determined by the cur- 
vature of the dispersion around the lowest excitation, \i 
the chemical potential given by n = g^B,{H — H c ), U the 
interaction constant and iV the number of dimers. If a 
parabolic isotropic dispersion relation = (fifc) 2 /2m is 
used, then the critical chemical potential is given by 
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This relation leads to the phase boundary described by 
the power law 

(g/2)[H N (T)-H c ] = AT*, (3) 

with critical exponent </>bec — 3/2, 2 where Hn(T) is 
the transition field at temperature T. A point given by 
T = and H = H c on the temperature vs field diagram 



denotes the quantum critical point (QCP). Equation (2) 
or (3) gives the critical behavior near the QCP character- 
istic of the magnon BEC. The BEC of magnons has been 
studied extensively in many gapped spin systems 5-12 and 
the power law behavior of the phase boundary was con- 
firmed. 

TICUCI3 is an S = 1/2 interacting spin dimer sys- 
tem in which a chemical dimer CU2CL3 forms a antifer- 
romagnetic spin dimer. The interactions between neigh- 
boring dimers are three-dimensional. The lowest excita- 
tion occurs at Q = (0,0,1) and its equivalent recipro- 
cal points. 13-15 The magnitude of the excitation gap 
is A/fce = 7.5 K. 4,16 In the previous magnetization 
and specific heat measurements in magnetic fields on 
TICUCI3 , 4, 17 it was shown that the phase boundary for 
the FIMO is expressed by the power law with critical ex- 
ponent (f> = 2.0 ~ 2.2. This critical exponent is somewhat 
larger than </>bec = 3/2. However, in a lower tempera- 
ture region 0.5K < T < 3K Shindo and Tanaka obtained 
rather smaller value <f> = 1.67 by the specific heat mea- 
surement in magnetic fields . 18 

Recently, the deviation of (f> toward larger value from 
0bec = 3/2 has been discussed theoretically . 19-22 Using 
stochastic series expansion quantum Monte Carlo sim- 
ulations, Nohadani et al. 19 studied the FIMO in cubic 
lattices of dimers with antiferromagnetic Heisenberg in- 
teractions between dimers. They showed that <p decreases 
with decreasing fitting temperature range and converges 
at 4>beg — 3/2, and that when fitting range is small 
enough, (f> is independent of details of the system. They 
ascribed this behavior to temperature-driven renormal- 
ization of the quasiparticle effective mass and of the ef- 
fective chemical potential. On the other hand, Sherman 
et al. 20 argued the critical exponent 4>, assuming the rel- 
ativistic dispersion of the form Ek = \/A 2 + Ck 2 , be- 
cause the real dispersion curve around the lowest exci- 
tation is better described by the relativistic form than 
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the parabolic one. They showed that the deviation of 
the dispersion curve from the parabolic form gives rise 
to the experimental larger exponent, and that the ex- 
ponent converges at </>bec = 3/2 with decreasing fitting 
range. Kawashima 21 demonstrated analytically and nu- 
merically that the critical exponent 0bec = 3/2 derived 
by the HF approximation is exact, although the HF re- 
sult is usually incorrect for the critical behavior. Mis- 
guich and Oshikawa 22 extended the HF calculation by 
Nikuni et al. , 2 using a realistic dispersion relation 13, 15 
and achieved remarkable quantitative agreement with 
the experimental phase diagram. They also calculated 
the critical density of magnons n cr , which corresponds 
to the absolute value of the magnetization at Tn . 

To the best of the authors' knowledge, TICUCI3 is the 
best system to study the magnon BEC and the critical 
behavior, because details of the magnetic excitations are 
known and the BEC state can be reached using conven- 
tional superconducting magnet. To obtain precise critical 
exponent <f> and critical density n cr , we carried out mag- 
netization measurements on TICuCLj down to 77 mK. 

Single crystals of TICUCI3 were grown by the vertical 
Bridgman method. The details of preparation were re- 
ported in Ref. 4. The magnetization measurements were 
performed using SQUID magnetometer (Quantum De- 
sign MPMS XL) in the temperature region 1.8K < T < 
100K in magnetic fields of up to 7 T. The magnetic fields 
were applied parallel to the 6-axis and [2, 0, 1] direction 
and perpendicular to the (1, 0, 2) plane. The magneti- 
zation measurements were also performed using Faraday 
Force Magnetometer 23 at Institute of Solid State Physics 
in the temperature region 77 mK < T < 4 K in mag- 
netic fields up to 7 T using a dilution refrigerator. The 
magnetic field was applied perpendicular to the (1, 0, 2) 
plane. The gradient field of 5 and 8 T/m was applied to 
produce the Faraday force. 

Figure 1 shows low-temperature magnetization M 
measured at various external fields for H ± (1,0,2). As 
temperature is decreased, the magnetization decreases 
with convex function of temperature and then increases 
exhibiting the cusplikc minimum indicative of the three- 
dimensional magnetic ordering. This magnetization be- 
havior was observed, irrespective of field direction. We 
assign the temperature giving magnetization minimum 
to the ordering temperature Tn- For T > Tn, the magne- 
tization corresponds to the number of thermally excited 
magnons, which decreases with decreasing temperature. 
The BEC of magnons occurs at T = Tn, and below Tn, 
the number of condensed magnons increases with lower- 
ing temperature. The increase of the condensed magnons 
surpasses the decrease of thermally excited magnons. For 
this reason, the magnetization has the minimum at Tn. 
With decreasing magnetic field, Tn decreases, and the 
cusplike anomaly due to the phase transition becomes 
smaller. This is because the number of magnons relevant 
to the BEC decreases when the magnetic field approaches 
the critical field H c . 

The critical density n cr corresponds to the abso- 
lute values of the magnetization at Tn and is obtained 
from the minimum of the magnetization M cr as ti C y — 
(2/gftB)M a . We corrected the magnetization for the 




Fig. 1. (Color online) The temperature dependence of the mag- 
netizations of TICUCI3 measured at various magnetic fields for 
flX(l,0,2). 




Fig. 2. (Color online) Temperature dependence of critical den- 
sity of magnons in TlCuCl 3 for H\\b, H||[2,Q, 1] and i2"J_(l, 0, 2). 
The solid line is the theoretical calculation by Misguich and Os- 
hikawa . 22 



paramagnetism due to impurities which produce the 
magnetization described by the Brillouin function of 
H/T, and also for the Van Vleck paramagnetism and the 
diamagnetism due to core electrons. The g-factors used 
are g = 2.06 for H \\ b and H || [2,0, 1] and g = 2.23 for 
H _L (1,0,2), which were determined by ESR measure- 
ments. Figure 2 shows the critical density func- 
tion of temperature. The solid line in Fig. 2 is the critical 
density calculated by Misguich and Oshikawa 22 using the 
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Fig. 3. (Color online) Transition field vs critical density n cr 
obtained for H || b, H || [2,0, 1] and H _L (1,0,2). Solid lines are 
fits by eq. (4). 



realistic dispersion. 13,15 The experimental and theoret- 
ical critical densities agree well in the low temperature 
region of T < 3 K. However, in the high temperature 
region, the experimental values are larger than the theo- 
retical ones. The disagreement increases with increasing 
temperature. 

In the HF approximation, the transition field H^(T) 
is expressed as 22 



(g/2)[H N (T) - H c ] = 2Un CT {T). 



(4) 



Using eq. (4) and the experimental critical density shown 
in Fig. 2, we evaluate the magnon-magnon interaction 
constant U. Figure 3 shows the plots of transition field 
Hn(T) vs critical density n cr (T), where Hjq(T) is normal- 
ized by the ^-factor. The linear relation between H-^(T) 
and n cr (T) holds for n cr (T) < 3 x 10~ 3 . Fitting eq. (4) in 
this small n CI (T) region, we obtain U/Icb = 312 K, 311 
K and 315 K for H \\b,H± (1,0,2) and H || [2,0, 1], re- 
spectively. Their average is U/kB — 313 K. This value of 
U is somewhat smaller than U /fee = 340 K obtained by 
Misguich and Oshikawa 22 using the magnetization data 
reported in Ref. 4. 

To determine the transition fields for T < 1.8 K, we 
performed magnetization measurements using the Fara- 
day Force Magnetometer. Fig. 4 shows the magneti- 
zation curve measured between 77 mK and 1.2 K for 
H _L (1, 0, 2). The raw magnetization has the small back- 
ground in the low field region due to impurity phase 
which should be hydrate on the sample surface. In Fig. 
4, the background magnetization was corrected. The 
magnetization is almost zero up to the transition field 
i?N — 4.7 T and increases rapidly. The gapless ordered 
state appears for H > H^. The magnetization does not 
show sharp bend at , but is rather rounded even at 77 
mK. This is ascribed not to the gradient field to produce 
the Faraday force, because the distribution of magnetic 
field in the sample is less than 0.04 T. We infer that the 
antisymmetric interaction of the Dzyaloshinsky-Moriya 
type which mixes the singlet and triplet states gives rise 




H(T) 

Fig. 4. (Color online) Magnetization curves in TICUCI3 measured 
by Faraday Force Magnetometer for H _L (1, 0, 2) at various tem- 
peratures. Backgrounds due to impurities were subtracted. Each 
plot is shifted vertically by 2 X 10 — 4 /ib for clarity. Arrows denote 
transition fields iir-j(T). 



to the smearing of the magnetization anomaly at H^- We 
assign the transition field (T) to the field of inflection 
in the derivative of magnetization dM/dH, as shown in 
the inset of Fig. 4. 

The phase transition points obtained by temperature 
and field scans of magnetization are summarized in Fig. 
5. Since the phase boundaries for H || b, H \\ [2, 0, 1] and 
H _L (1,0,2) coincide when normalized by the g-factor, 
we can deduce that the phase boundary is independent of 
the external field direction, which implies that the mag- 
netic anisotropy is negligible in TICUCI3. Solid line in 
Fig. 5 is the HF calculation using the realistic disper- 
sion 13,15,22 and the interaction constant of U/k-B — 313 
K obtained by our magnetization measurement. The ex- 
perimental phase boundary for T < 5 K agrees almost 
perfectly with the calculation. 

We analyze the phase boundary for the magnetic field 
perpendicular to the (1,0,2) plane using the power law 
given by eq. (3). We fit eq. (3) in the temperature range 
of T m i n < T < T max , setting the lowest temperature at 
Tmin — 77 mK and varying the highest temperature T max 
from 4.6 K to 1.94 K. The critical exponent </> obtained 
by the fit is graphed in Fig. 6. It is observed that the crit- 
ical exponent 4> decreases with decreasing T max , and con- 
verges at 0bec = 3/2, as recent theory predicts. 19,20,22 
For T max < 2.4 K, we obtain (g/2)H c = 5.27 ± 0.01 
T, A = 0.225 ± 0.01 T/K^ and cj> = 1.50 ± 0.06. When 
(p = 3/2, the parameter A corresponds to the coefficient 
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Fig. 5. (Color online) Phase diagram normalized by the g-factors 
in TICUCI3. The solid line is the result of calculation using the 
realistic dispersion and the interaction constant of U/ks = 313 
K. 22 




Fig. 6. (Color online) Critical exponent </, as a function of T max 
obtained by fitting eq. (3) to transition points between 77 mK 
and T max (K). The transition points were collected for H _L 
(1,0,2) 



of T 3 / 2 in eq. (2), and is given by 
2.61 U fmkuX 



A = 



(5) 



Substituting A = 0.225 T/K 15 into eq. (5), we ob- 
tain mk^jf? = 0.0204 K _1 . This value is consistent 
with mfce/ft 2 = 0.0229 K _1 evaluated from the curva- 
ture of the dispersion relation around the lowest exci- 
tation. 22 In eq. (2), the lattice constant is chosen to be 
unity. If we use an average lattice constant a given by 
a = {abc sin /3) 1 / 3 , where a = 3.98 A, b = 14.14 A, 
c = 8.89 Aand (3 = 96.32° are lattice parameters in 
TlCuCl 3 , 24 we obtain m = 2.61 x 10~ 26 g. This effec- 
tive magnon mass is approximately 1/60 of the proton 



mass. This smaller magnon mass enables the BEC at he- 
lium temperatures in spite of small density of order of 
10" 3 . 

In conclusion, we have presented the results of mag- 
netization measurements performed on TICUCI3 at tem- 
peratures down to 77 mK and in magnetic fields up to 
7 T for three different field directions. The critical den- 
sity of magnons as a function of temperature and the 
magnetic phase diagram for external field vs temper- 
ature were obtained. The magnon-magnon interaction 
constant was estimated as U /ks = 313 K. The phase 
boundary is expressed by the power law of eq. (3) and 
agrees almost perfectly with the magnon BEC theory 
based on the Hartree-Fock approximation with realistic 
dispersion relations 22 and U/ks = 313 K obtained by 
our magnetization measurement. The critical exponent <f> 
decreases with decreasing fitting temperature range. For 
77 mK < T < 2.4 K, we obtained (j> = 1.50 ± 0.06, which 
coincides with 0bec = 3/2 derived from the magnon 
BEC theory. 2 ' 21 The effective magnon mass obtained by 
the present analysis is consistent with that evaluated 
from the curvature of the dispersion relation around the 
lowest excitation. The present results strongly support 
the BEC description of the field-induced magnetic or- 
dering in TlCuCl 3 . 
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